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ABSTRACT 

A first adiabatic core is a transient object formed in the early phase of star 
formation. The observation of a first core is believed to be difficult because of its 
short lifetime and low luminosity. On the basis of radiation hydrodynamic simu- 
lations, we propose a novel theoretical model of first cores, Exposed Long-lifetime 
First core (ELF). In the very low-mass molecular core, the first core evolves slowly 
and lives longer than 10,000 years because the accretion rate is considerably low. 
The evolution of ELFs is different from that of ordinary first cores because radi- 
ation cooling has a significant effect there. We also carry out radiation transfer 
calculation of dust- continuum emission from ELFs to predict their observational 
properties. ELFs have slightly fainter but similar SEDs to ordinary first cores 
in radio wavelengths, therefore they can be observed. Although the probabilities 
that such low mass cores become gravitationally unstable and start to collapse 
are low, we still can expect that a considerable number of ELFs can be formed 
because there are many low-mass molecular cloud cores in star-forming regions 
that can be progenitors of ELFs. 

Subject headings: stars: formation — ISM: clouds — radiative transfer — hy- 
drodynamics 



Introduction 



Larson! (119691 ) studied star-formation processes in low-mass molecular cloud cores us- 



ing one-dimensional (ID) spherically symmetric hydrodynamic simulations. In his work, 
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a transient, quasi-hydrostatic object is formed in the early phase of the proto-stellar col- 
lapse, a so-called first (adiabatic) core. The first core evolves via accretion from the natal 
core and it collapses again when it attains a temperature high enough to dissociate hydrogen 
molecules. Although first cores are transient, they are essential to understand star formation 
processes because they assume a crucial role in the two major problems in the star formation 
process, the angular momentum and the magnetic flux problem. The first cores are inter- 



esting sites of diverse phenomena related to t 


lese problems such as driving molecular out- 


flow 


Tomisaka 20021; Matsumoto & Tomisaka 


2004: Machida et al. 20061: Duffin & Pudritz 


2009; 


Tomida et a . 


2010l 


Ciardi k, Hennebelle! 


20101). binarv formation (Machida et al. 


20081; 


Commercon et al. 


2010). planet formation (Inutsuka et al. 2010). circumstellar disk forma- 


tion 


(Hennebelle & Ciardi 2009: Machida et al. 2010: Machida & Matsumoto 2010: 


Bate 


2010; 


Dapp & Basi 


l 2010 


), and so on. More than 40 years have passed since Larson's the- 



oretical prediction, but it has not been robustly confirmed yet by observations due to its 
observational difficulties. The lifetime of a first core is so short, on the order of 100-1000 
yrs, compared to the dynamical time-scale of its natal cloud, which can be about 10 4 — 10 6 
yrs depending on its initial density, that the detection probability of the first core is consid- 
ered to be small. The first core is also deeply obscured in the infalling envelope with very 
large column density; hence it is believed to be difficult to observe with current instruments. 
However, if we could observe first cores directly, they will significantly contribute to our 
understanding of star formation processes. 



The evolution of cloud cores with ma sses around 
dimensional (3D) numerical simulations (Bate! 19981 



20081 : ICommercon et al.lhoio l: iTomida et aDboiol : featel boio L In those cases, the first cores 



M(p has been well studied using three- 



Whitehouse fc Bate! 120061 : ISaigo et al. 



evolve over very short dynamical time-scales because accretion from their natal cores con- 
trols the evolution. Rotational support extends the lifetim e of the first cores , but it is still 
brief compared to the dynamical time-scales of the clouds ( ISaigo et al.l 120081 ) . However, we 
can expect that the first cores formed in very low-mass cloud cores do not evolve into the 
second collapse stage over these short dynamical time-scales because most of the gas in the 
envelope quickly falls into the first core and the accretion rate rapidly decreases. These first 
cores survive longer before the second collapse starts than do cores formed in more massive 
clouds. In the long-term evolution of these systems, ra diation cooling plays an important 
role in disk stability and angular momentum transport ( iGammid 120011 ) . To handle the gas 
thermodynamics properly, radiation hydrodynamics (RHD) simulations are required. 



Very recently, several groups reported possible detection of first core candidates. IChen et al 



( 120101 ) found a very faint and compact core and claimed that it was a good first core can- 
didate. Their object must be very young, but it seems to be more evolved than a first core 
because it is associated with a well-collimated and high-velocity outflow, which is thought to 



-3 - 



be driven from a protostar after the second collapse. IChen fc Arcd (120101 ) reported another 
candidates in the R Corona Australis cloud. Enoch et al.l (120101 ) found a good candidate for a 
first core whose spectral energy distribution (SED) is consistent with a theoretical model. A 
Japanese group (Kawabe et al., in preparation) also reported a peculiar dense core which is a 
good first core candidate. These objects need to be confirmed by more precise observations, 
but these observations imply that first cores are more common than expected. The observed 
candidates commonly have faint SEDs and some of them seem to have small masses. Con- 
sidering observed c ore mass function ( CMF), there are so many low-mass cores whose masses 



are le ss than 1M (IMotte et al.lll998l : Enoch et al.ll2007l : iRathborne et al.ll2009l : lAndre et al. 



2010). In the initial m ass function (IMF), there are still many low-mass stars and brown 
dwarfs ( jChabrierl 120031 ) . although the relationship between the two mass functions is still 
unclear. Therefore it seems reasonable and interesting to study the evolution of low-mass 
cores which have masses less than 1M . 

The observational p roperties of first cores such as their SEDs were predicted on the basis 
of ID RHD simulations (IMasunaga fc Inutsukall2000l ) and post-processing radiation transfer 
calculati ons by using the results of 3D hydrodynamic simulations with barotropic approx- 
imation (ISaigo fc Tomisakall2010l ). First cores have faint, low-temperature black-body-like 
SEDs in wide wavelengths from mid-infrared to centi-meter-length radio wave. First cores 
are faint but detectable with modern facilities such as the Herschel Space Observatory, 
the Submillimeter Array (SMA), the Expanded Very Large Array and the Atacama Large 
Millimetre/sub-millimetre Array (ALMA). With ALMA in full operation, it is expected that 
we will be able to resolve the structure of first cores such as spiral-arms formed via gravita- 
tional instability in the disk and driving region of the molecular outflow driven by magnetic 
fields predicted with magnetohydrodynamics (MHD) simulations. 

In this paper, we propose a new scenario for the evolution of first cores formed in 
very low-mass cloud cores. We prove that first cores in the very low mass molecular cloud 
cores have significantly long lifetimes. We also calculated SEDs and visibility amplitude 
distributions of our model, and compared them with the models of both a typical first core 
and a starless core. The plan of this paper is as follows. In section 2 we briefly describe our 
method of RHD simulations and models used in the simulations. We show the results of the 
numerical simulations in section 3 and the observational properties of those simulated first 
cores in section 4. Section 5 is devoted to conclusions and discussions. 



-4- 



Method and Models 



We perform 3D nested-grid self-gravitational RHD simulati ons. Radiation transfer 



is tre ated with gray flux limited diffusion (FLD) approximation (jLevermore fc Pomraning 



1981). The details of our simulation code are described in our previous work ( ITomida et al 



20101 ) . but we neglect magnetic fields in this work because initially we just intend to prove 
the concept. We use an idealized equation-of-state (EOS) with adiabatic index 7 = 5/3 
which is different from that used in our previous work, 7 = 7/5, to make t he calculations 



easier . We use the compiled table s of Rosseland and Planck mean opacities of ISemenov et al. 



( 120031 ) and iFerguson et al.l ( 120051 ) . as in our previous paper. The number of grid points in 
each level of the nested-grids is 64 3 , and the finer grid is generated around the center of the 
computational domain to resolve the local Jea ns length with at least 16 grid cells to prevent 
artificial fragmentations ( ITruelove et al.l 1 19971 ). 



the center and the radii of the clouds are (p c ,R c ) = (3.2 x 10 18 gem 3 
and Rl, and (p c ,R c ) = (3.2 x 10~ 16 gem" 3 



We calculated three models for comparison: SI is a non-rotating 1M model, Rl is a 
rotating 1M model, and R01 is a rotating very low- mass model whose mass is O.1M . As 
the initial conditions, we take critical Bonnor-Ebert-like spheres with uniform rotation, and 
increase the density by a factor of 1.6 to make them unstable. The initial gas densities at 

6300 AU) in SI 
630 AU) in R01. It is not trivial how we should 
scale the angular velocity between models with different masses, but here, we assume that 
both Rl and R01 initially have the same amount of rotational energy, T = aMR 2 Q 2 (a is a 
constant of order unity). The initial angular velocities are 4.3 x 10~ 14 s _1 and 1.4 x 10~ 12 s~ x 
in Rl and R01, respectively. We adopt the boundary conditions that all the cells outside 
the sphere of the cloud radius maintain their initial values. The initial gas temperature and 
boundary conditions for radiation transfer are set to 10K. Here, we note that our boundary 
conditions allow the gas to inflow into the computational domain through the boundaries, 
and at the end of the simulations about 30% of the initial mass increased in the low-mass 
case and less than 10% increased in the 1M Q cases. 



3. Results 

We show the cross sections of the gas density and temperature of Model Rl and Model 
R01 at the epoch when they have the same first core masses, 1.06 x 10 _1 M Q , in Figure [TJ 
Here, the first core mass Mpc is defined as 



M FC = P AVc > 

P>P¥C 



(1) 
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where AV C is the volume of the cell and pfc is the critical density defined as the minimum 
density in the gas that experienced the shock, which we identified as the region where the 
radial velocity is smaller than the local sonic speed, \v r \ < c s . The age of the first core at 
this epoch, t-pc, is 3,100 yrs in Rl and 10,600 yrs in R01 from the core-formation. Model 
R01 has a significantly larger (~ 100 AU) first core disk because it is more evolved and the 
gas with a larger specific angular momentum from the outer region has accreted onto the 
core. Another outstanding difference between the models can be seen in the gas density and 
temperature in their envelopes. For the low mass core model, almost all the materials in 
the natal core have already accreted onto the first core, and therefore, the gas density in the 
envelope is reduced drastically. Two factors cause the lower temperature in the envelope of 
R01: the smaller accretion rate in the very low-mass core results in the weaker shock at the 
surface of the first core, and thus, less entropy is produced at the shock (in other words, the 
first core is intrinsically colder and fainter), and the smaller optical depth of the envelope 
contributes to the efficient radiation transport and cooling. 

In Figure [2j we show the time-evolution of physical quantities in each model such as 
central gas density (a), temperature (b), first core mass (c) and accretion rate onto the 
first core (d). It is obvious that R01 has a smaller accretion rate and its mass increases 
more slowly than it does in other models. As a result, the gas density and temperature at 
the center of the first core evolve significantly slower. The lifetime of the first core, from 
its formation to the second collapse, reaches more than 10 4 yrs in R01. Both 1M© models 
achieved similar accretion rates because the initial structure of the cloud core determines the 



accretion rate (ISaigo et al.ll2008[ ). although the non- rotating model evolves faster than the 
rotating model. These results clearly indicate that centrifugal force supports a considerable 
mass and prevents the first core from colla psing. We note that the lifetim e of the non-rotating 



model is longer than previous predictions (IMasunaga fc Inutsukall2000l ). by about a factor of 



two, because we assume a simple stiff EOS of 7 = 5/3, and our model has a smaller accretion 
rate due to stable initial conditions. But in the rotating models, the difference of EOS affects 
the results less significantly because the centrifugal force dominates the structure of the first 
core disk. 

We show the thermal evolution tracks of the gas elements at the center of the clouds 
in the p — T plane in Figure |3J All the models show the effects of radiation cooling, but it 
appears most significantly in the low mass model. We also plot the evaporation temperatures 
of each dust component, which affect the thermal evolution of the gas. When the gas 
temperature reaches the evaporation temperature of iron and silicates, T ~ 1400 K, the 
opacity drops substantially and the core collapses violently, similar to the second collapse 
due to the endothermic reaction of hydrogen molecule dissociation. However, if we adopt 
a soft EOS with 7 = 7/5, the impact of dust evaporation becomes less important because 
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the gas density (and therefore the optical depth) at the same temperature is higher. Thus 
the effect of radiation cooling is important in the low mass core where dynamical accretion 
does not dominate the evolution. RHD simulations are required for these systems because 
the barotropic approximation fails to reproduce the realistic thermal evolution. 

Our results show that the low-mass model has a longer lifetime compared to the dy- 
namical time-scale of its natal core. This suggests that a large fraction of low-mass cores 
can harbor the first core when we carry out an unbiased survey of collapsing starless cores. 
Although only a small fraction of the low-mass cloud cores may be gravitationally bound 
and will collapse into stars, there are a lot of low mass cores in star-forming regions, and we 
can, therefore, expect that a considerable number of ELFs can form. 



4. Observational Properties of the First Cores 

We calculated the SEDs of first cores in the RHD simulations by performing post- 
processing radiation transfer. We solve the following radiation-transfer equation in one di- 
rection to draw an intensity distribution map for various inclination angles and wavelengths: 

^= P k v (B v (T)-I v ), (2) 

where B V {T) is the Planck function of temperature T and k v is the monochromatic (absorp- 
tion) opacity. Here, we assume that the temperature T obtained in the RHD simulations 
is correct, and we ignore scattering, which is a valid assumption in mid-infrared or longe r 



wavelengths. We adopt the monochromatic dust opacities provided by lSemenov et al.l (120031 ). 
using the homogeneous aggregates model of normal abundances. 

In Figure HI we show the SEDs of our RHD models in the face-on configurations. For 
comparison, we also plot the SED of O.1M Bon nor-Ebert sphere, a model of a low mass 



starless core, and observed SED of L1521F-IRS (IBourke et al.ll2006l ). which is a very low- 
luminosity object (VeLLO) in the Taurus molecular cloud. The distance towards the targets 
is set to 150 pc and the SEDs are measured with a (1000 AU) 2 aperture. The flux in the far- 
infrared wavelengths increases when the first core forms and it can be observed with Herschel. 
Compared to L1521F-IRS, Model Rl has similar brightness in submillimeter region but it 
disagrees in mid-infrared region. In contrast, Model R01 is fainter than L1521F-IRS in all 
the wavelengths. In radio wavelengths, in spite of having a mass 10 times larger than that 
of R01, Rl is only about twice as bright as R01. However, the peak flux around the far- 
infrared wavelengths is significantly larger in Rl. This is because it has a warmer first-core 
disk and envelope, as we mentioned before (Fig. [[]). In contrast, R01 is more luminous in 
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the mid-infrared wavelengths because it has a thinner envelope. If we follow the more long- 
term evolution of the low mass cloud model, this tendency is enhanced. Our models may 
explain the weak em ission in mid-infrared wavelengths in some observed first core candidates 



( [Enoch et al.ll2010l ). We will discuss the time-evolution of observational properties of first 



cores in a subsequent paper. 

First cores emit larger flux in the infrared region than starless cores, but it is still 
difficult to identify the existence of a first core only from the SED, because it is so faint 
in the infrared wavelengths and the differences in radio wavelengths between the first cores 
and the starless cores are not so prominent. Another good method to distinguish first cores 
from starless cores is high resolution observation measuring visibility amplitude (Fourier 
components of the intensity map) distributions with (sub)millimeter interferometers. We 
show simulated visibility amplitude profiles of Model R01 measured in 850 /im, obtained 
with the Common Astronomy Software Applications (CASA) in Figure [5j The visibility 
amplitude of the starless core decreases steeply in the small scale because it contains no fine 
structure. In contrast, the first core clearly shows a shallow distribution. This feature is 
not solely seen in the ELF but common in first cores. This is firm evidence of the existence 
of the first cores. Current observations with SMA is not sufficient to resolve the first core 
directly, but sufficient to identify the first core. ALMA in full operation will reveal the 
detailed structure of the first core. 



5. Conclusions and Discussions 

We performed 3D RHD simulations of low-mass cloud cores and showed that first cores 
formed in very low mass (O.1M ) cloud cores live more than 10 4 yrs. Those first cores 
have thin envelopes, therefore we name this first core an "Exposed Long-lifetime First core 
(ELF)." ELFs experience different evolution from ordinary first cores whose evolution is 
dominated by accretion from the natal core. ELFs use up the gas in the envelope in the 
early stages of its evolution, and then the mass and angular momentum redistribution in 
the disk control the evolution of ELFs. Radiation cooling plays a critical role there similar 
to the Kelvin-Helmholtz contraction, dominating the disk stability and angular momentum 
transport with the time-scale longer than the dynamical time-scale of accretion. Thus, the 
evolution of ELFs is different from that of the first cores formed in ordinary mass cloud 
cores, which have been well studied so far. The gas in ELFs does not behave adiabatically 
anymore and the barotropic approximation breaks down. 

We also calculated the observational properties of ELFs. We found that in radio wave- 
lengths (300 fim ~) they have slightly fainter but quite similar SEDs to ordinary mass first 
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cores. This fact suggests that ELFs are detectable as well as ordinary first cores and they 
can be detected even with current instruments like SMA. On the other hand, ELFs are 
more luminous in mid-infrared region because they have less massive envelopes than those 
of typical first cores. 

Very low-mass cores will not collapse so often by self-gravity and the formation proba- 
bility of ELFs may be low. However, ELFs are stil l worth cons i derin g; because there are 



a number of low-mass cores both in observations ( Motte et al.l 11998 



Rathborne et al. 



2002; 



20091: 



Hennebelle fc Chabrier 



Andre et al.ll2010l) and in theoretical predictions ( Padoan fc Nordlund 



Enoch et al. 



2007; 



20081 ). ELFs are long-lived compared to the dynamical time- 



scale of their natal cloud cores, so observation possibilities are higher than for usual first 
cores. Therefore, we can expect a considerable number of ELFs exist in the star-forming 
regions. Ot her mechanisms such a s external pressure, cloud-c loud collisi on, radiation driven 
implosion (IMotoyama et al.l 120071 ). and dynamical ejection (iBatd 120091 ) may help the for- 
mation of ELFs. We will study the effects of magnetic fields and realistic EOS in future 
works. 
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Fig. 1. — Horizontal (a, d) and vertical (b, c, e, f) cross sections of gas density (a, b, d, e) and 
temperature (c, f) at the epoch of the same first-core masses, (a)-(c): Rl at t FC = 3, 100 yrs, 
(d)-(f): R01 at tpc — 10,600 yrs. The first cores are indicated by white dashed lines. 
Pfc = 1.0 x 10 -14 g cm -3 and 7.6 x 10~ 16 g cm -3 in Rl and R01, respectively. Despite the 
same first core masses, the envelope in R01 is clearly depleted and colder. Note that the 
spatial scales in (a) and (b) are different from others. 



-12- 




Time [yr] 



Fig. 2. — Time-evolution of the physical quantities, central gas density (a), temperature (b), 
first core mass (c), and smoothed accretion rate (d). Model R01 shows prominently longer 
lifetime than 1M Q models, SI and Rl. 
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Fig. 3. — The evolution tracks of the thermal properties of the central gas elements in the 
p — T plane. The dust evaporation temperatures are over-plotted with orange dash-dotted 
lines. All the models show the influence of radiation cooling but it is most significant in the 
low mass model. 
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Fig. 4. — SEDs of the first core models at the same epoch as in Figure Q] i n face-on config- 
uration and a 0.1M Q star-less core. The observed SED of L1521F-IRS (IBourke et al.ll2006l ) 
is also plotted. The distance towards the targets is 150 pc and the aperture is (1000 AU) 2 . 
First cores are more luminous than starless cores, especially in the far-infrared wavelengths. 
Model Rl is brighter than Model R01 in radio wavelengths. On the other hand, Model R01 
exceeds Model Rl in the mid-infrared region. 
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Fig. 5. — Visibility amplitude distributions of the 0.1 M first core model and a 0.1 M 
Bonnor-Ebert sphere as a model of a starless core. The first core model shows clearly 
shallower distribution compared to the starless core. The edge-on configuration shows more 
widely scattered visibility amplititude corresponding to its oblate morphology. 



